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ABSTRACT 


A transient ballistics and combustion model is derived to represent the 
closed vessel experiment that is vndely used to characterize propellants. The 
model incorporates the nitramine combustion mechanisms which are contained within 
the steady-state model developed as a part of last year's (FY 1976) research. A 
computer program is developed to solve the time-dependent equations, and is 
applied to explain aspects of closed vessel behavior. It is found that the rate 
of pressurization in the closed vessel is insufficient at pressures of interest 
to augment the burning rate by time-dependent processes. In the case of nitramine 
propellants, however, the cratering of the burning surface associated with combustion 
above break-point pressures augments the effective burning rate as deduced from 
the closed vessel experiment. Low pressure combustion is significantly affected 
by the ignition process and, in the case of nitramine propellants, by the developing 
and changing surface structure. Thus, burning rates deduced from the closed vessel 
experiment may or may not agree with those measured in the equilibrium strand 
burner. Application of the model to closed vessel test cases shows good agreement 
between theory and experiment. 

Series of T-burner experiments are performed to compare the combustion instability 
characteristics of nitramine (HMX) - containing propellants and ammonium perchlorate 
(AP) propellants. It is found that the inclusion of HMX consistently renders the 
propellant more stable. Although ash produced by more fuel-rich propellants could 
have provided mechanical suppresssion, results from clean-burning propellants permit 
the conclusion that HMX reduces the acoustic driving. 

Additional strand burning rate data are presented which provide further confirma- 
tion of the results of the steady-state modeling effect o^ FY 1976. 
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SECTION 1 


OBJECTIVE 

The objective of this research program is to study two aspects of the transient 
combustion properties of nitramine smokeless propellants. One is to determine the 
role of the nitramine ingredient in the combustion driving of acoustic instability 
in solid rocket motors. The second is to determine the effect of a transient 
pressure environment on the combustion of nitramine propellants in armament or 
closed vessel apparatus. The steady-state combustion characteristics of nitramine 
propellants have been described as a result of work performed under AFOSR Support 
Agreements AFOSR-ISSA-75-0005 and AFOSR-ISSA-76-0006 (Ref. 1). The present research 
begins to investigate unsteady-state or transient combustion characteristics. 


SECTION 2 


INTRODUCTION 

Nitramines are of interest for gun and rocket propellant applications because 
they are an energetic source of smokeless combustion products (Refs. 2-4). Recent 
studies have been devoted to understanding and improving the steady-state combustion 
characteristics of nitramine propellants (Refs. 1, 5-10). Rocket motors, are, 
generally, principally concerned vn'th steady-state performance. Although it is 
recognized that gun propellants perform in a highly transient pressure environment, 
there is basis for the inference that improvements in steady-state combustion will 
be reflected as improvements in the transient ballistics performance (Ref. 11). 

Anomalous combustion behavior can be brought about by unsteady processes occur- 
ring in combustion chambers. In rocket motors, the most prevalent behavior is com- 
bustion instability. Combustion instability is of particular concern in smokeless 
propellants because the absence of aluminum pov/der deprives the system of significant 
stabilization (Ref. 12). Virtually no work as been done to describe the influence of 
nitramines on combustion instability properties. In projectiles, the extreme trans- 
ient pressure environment can alter the combustion process from the steady-state as 
a result of time lags in the combustion response. Although this problem has been 
studied extensively with respect to rocket motor ignition and controlled termina- 
tion (e.g., Refs. 13-17), it has not been applied to gun propellants or to the more 
extreme projectile environment. Burning rates derived from closed-vessel (transient 
pressure) testing differ from those measured at constant pressure, which raises a 
question about the predictability of ballistics performance (Ref. 18). Further, the 
transient pressure environmenb may (Ref. 11) or may not (Ref. 19) aggravate any 
steady-state combustion anomaly to impair the quality of the impetus delivered. 

It is considered that understanding the transient combustion characteristics of 
nitramine smokeless propellants will provide additional direction for future develop- 
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mental efforts. The mechanistic understanding dervied from the steady-state research 
provides an essential background, and an understanding of the transient phenomena 
vnll afford improved tailoring and predictive capabilities. 


SECTION 3 


TRANSIENT COMBUSTION MODEL DEVELOPMENT 

3.1 REVIEW OF CLOSED VESSEL BURNING RATE DATA 

A reviev/ was conducted of available data wherein direct comparison mre mads 
between burning rates as measured in a closed vessel (transient pressure) and as 
measured in a strand burner (equilibrium pressure) (Refs. 20-24). Som.e of the closed 
vessel data v;ere accompanied by the associated test pressure-time, or pressure-dp/dt 
(P-dot), histories. As will be discussed below, it is important to know the P-dot 
values for transient analysis. Some of the data were for nitramine propellants; 
some were for homogeneous (unfilled single-base or double-base) active binder pro- 
pellants. In general, there appears to be very good agreement between the closed 
vessel and strand data for the homogeneous propellants. Such was the conclusion of 
Grollman (Ref. 23) in his study. Figures contained in Ref. (24) show a virtual 
identity between closed vessel and strand data. Lenschitz' data for homogeneous 
propellants (Ref. 20), some of which are reproduced here as Fig. 1, also show good 
agreement between the two methods. On the other hand, there is a significant and 
consistent difference betvjeen closed vessel and strand burning rates in the case of 
propellants containing HMX (Refs. 20-22). Examples from Ref. (20) are presented 
here as Fig. 2. At low pressure, the closed vessel rate tends to lag the strand 
rate; at high pressure, the closed vessel rate tends to lead the strand rate. As 
a result, the closed vessel data exhibit a more uniform and higher apparent pressure 
exponent. 

3.2 APPLICATION OF SIMPLIFIED "P-DOT'' THEORY 

VIhere there are differences between closed vessel and strand data, such 
differences have been attributed to the deductive nature of the closed vessel 
experiment. The closed vessel experiment measures a variable pressure versus 
time vjhereas the strand experiment measures directly the distance burned versus 
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time at a constant pressure. From the pressure versus time, ballistics analysis 
yields mass generation rate versus time If certain assumptions are made. If the 
geometry of the sample is known, the burning rate versus time may then be calculated. 
Thus, burning rate may be cross-plotted versus pressure. However, It has been 
suggested that the assumptions may be Incorrect and that the geometry of the sample 
may not be knov/n. Unknown factors frequently mentioned Include heat losses, non- 
uniform Ignition and burning, and charge break-up (Ref. 18). Low pressure uncertainties 
have limited the reporting of data to pressures above 5000 psi. 

A mechanistic basis for a difference between burning rates measured under 
steady-state and transient pressure conditions was developed in the course of 
studies of ignition and depressurization transients in solid rocket motors (Ref. 13). 
This mechanism was the thermal wave relaxation time In response to a pressure transient. 
A simplified form of expression was as follows (Ref. 14): 


r. 


1 = 1 + 2n^. P 

- 2 P 


H) 


r- = Instantaneous burn rate at pressure, P 

r = steady-state burn rate at P 
n = steady-state pressure exponent at P 
K = propellant thermal diffusivlty 
P = Instantaneous rate of change of pressure at P 

More detailed analysis resulted in the development of computer programs such 
as described in Ref. (16). Some of these programs have been reviewed by Kooker 
(Ref. 25) in the context of gun ballistics. 


PAtSE BLANK NC3T 
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An 'brder-of-magY)itude analysis of Eq. (1) reveals that P-dot effects should 
be small in the usual closed vessel experiment. Typical values of the parameters 
are as follows: 


r = 2 in. /sec @ 10,000 psi (Fig. 1) 
n = 0.8 (Fig. 1) 

p 

K =0.0002 in. /sec. (typical) 

P = 1000 P (Ref. 20) 

The impact of P-dot upon burning rate falls below 10% at pressures above 
10 Kpsi , and decreases with increasing pressures and burning rates. This result 
is consistent with the data in Fig, 1. Although the values of P-dot are very 
high when compared to values encountered in rocket motors, the values of pres- 
sure and r-squared also become very high so as to diminish the effect. Since 
most of the reported closed vessel data are at such high pressures, it would 
appear that agreement between the closed vessel and strand data should be 
expected assuming that P-dot would be the only reason for any difference. This 
appears to be borne out by the homogeneous propellants. 

Results of P-dot calculations applied to the nitramine propellant data are 
included as the dashed lines in Fig. 2. It v;as possible to do so because the 
P-dot histories for these tests v/ere available (Ref, 20). Note that the calcu- 
lated effect has the tendency to reduce the pressure exponent in the high slope 
region below 5 Kpsi. The calculated effect becomes small at high pressure, as 
discussed previously. However, the data behave quite differently from the cal- 
culated effect. In view of the fact that experimental uncertainties and P-dot 
effects appear to be adequately resolved in the case of homogeneous propellants, 
it follows that something different must be happening in the case of nitramine 
propellants. 
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3.3 THE ANALYTICAL M'^DEL 

3.3.1 Conceptual Basis 

The analytical model for the transient combustion of nitramine propellants 
is based upon the nitramine surface structure mechanisms contained v/ithin the 
steady -state model (Refs. 1, 7), coupled to a thermal wave relaxation model 
(Ref. 16) of a homogeneous solid. 

At low pressure, and following ignition, a planar melt layer must form 
and develop. This development will be an endothermic contribution to the 
energy balance, will tend to oppose the ordinary collapsing of the thermal wave 
during a pressure rise, and will tend to delay a break point. These effects 
should produce a burning rate lag at low pressure. At high pressure, and where 
coarse particles are present, the development of the surface craters in accordance 
with break points will increase the burning surface area and thereby augment the 
mass burning rate. As discussed earlier, the closed vessel deduces the mass 
burning rate and does not measure directly the linear burning rate. Therefore, 
if the burning surface is underestimated, the linear rate as computed from the 
closed vessel data will appear to lead. 

Examples of cratered surfaces appear in Fig. 9 of Ref. (1). If the cratering 
mechanism increases both the linear burning rate and the surface area of the 
small pellets used in transient devices, then the mass burning rate would be 
proportional to the square of this process. Thus, the exponent break point 
would have a more serious effect in pressurizing chambers. 

3.3.2 Basic Equations for a Homogeneous Propellant 

Instantaneous burning rate is related to a reference or steady-state value 
by means of an Arrhenius equation for surface decomposition: 
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(2) 




= r exp [(-E^/f^„)0/T_ - l/f_) ] 

s 0 _ JJ s ^ ^ 

Eg = activation energy of surface decomposition 
Rjj = universal gas constant 
Tg " instantaneous surface temperature 
Tg = surface temperature corresponding to r 
r = steady-state burning rate at P, known from data or 
the steady-state model 


The instantaneous surface temperature is calculated from the time-dependent 
Fourier equation: 


aT 
a t 



+ r. 




exp (-E 3 /RJ) 


(3) 


T = temperature distribution 
t = time 

X = distance into the solid 
Qgg = condensed phased heat release 
A = prefactor 

c = propellant heat capacity 


For numerical solution, Eq. (3) is written in a suitable finite difference 
form such that the nominal mesh size is selectable depending upon the thermal 
wave thickness. The selection affects accuracy, stability and computer time. 
The initial and boundary nditions for Eq. (3) are: 


OEIGWM' page & 
01? POOE quality 
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( 4 ) 




T (X.O) = T 


3T 


8X 


(w,t) = 0 


M- 


" I ? = W * »'■ <='^3 


To! 


«ss = *!s 


w = web thickness 

Tq = initial bulk temperature 

f t ■ external heat source, for example as would be 
input to produce ignition 

k = propellant thermal conductivity 
« 

Qg = heat flux to the propellant surface 
p = propellant density 

The pressure variations are computed by a mass and energy balance for the 
chamber: 

(V-nm^) ^ = (OTg + oP) -gr ’""c al 


Tg = gas temperature 
m^ = mass of gas stored in chamber 

R = gas constant 

V = chamber gas volume 

n = covolume 

dm .y 

ar ’ '’1 3t ' 

Aj^ = burning surface area 

For a closed vessel, it is assumed that there is no mass exhaust. 


( 6 ) 
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(7) 


dT 




-4t ‘ " '^g' 


dm. 


dt 



c 


g 


Y = ratio of specific heats 

Tp = adiabatic flame temperature {a function of P) 

c = gas heat capacity (a function of P) 

* 

= heat loss rate to wall 

V\ = "c'l'g - * (-)g 'V - C) (8) 

A, = wall area 
w 

= suitable heat transfer coefficient 
= wall temperature 

(ea)g = product of gas em1ss1v1ty and Stefan - Boltzmann constant 


The unsteady heat transfer in the wall Is solved In analogous fashion to 
Eqs. (3) and (4), using Eq. (8) to describe the flux at the boundary; but r^.=0 and 
Qgs~0 for the wall. 

The foregoing equations, plus Eq. (3) written for the wall, constitute a system 
for the unknowns P, V, m^, Tg, Q^, T^, Q^, T^ and r^. as functions of time. These 
equations are solved numerically by a forward difference method. The instantaneous 
burning rate will differ from the steady-state value at a given pressure to the extent 
that T^ differs from T . This difference can arise as a result of the changing thermal 
profile in the solid in response to the transient heating. 

3.3.3 Equations Related to the Surface Structure of Nitramine Propellants 

Additional equations enter into the system for nitramine propellants because 
of the relationship of the surface structure to the burning rate (Refs. 1, 7) and 
mass flow rate. It is assumed that the surface structure will develop in three 
stages in the course of the pressurization transient: 
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a) formation of the melt layer following ignition; 

b) disappearance of the melt layer when the break point criterion is 
reached; 

c) formation of the surface craters following disappearance of the melt 
layer. 

The melt layer will start to form as soon at Tg = the melting point of 
the nitramine. The significance of the melt layer thickness is that it positions 

the heat of fusion, which subtracts from the exothermic in the energy 
balance equations. It is expected that the location of T^^ will move into the solid 
with time; Tj^j is 551 °K for HMX, and computed steady-state surface temperatures are 
circa 1000°K. 

The second stage commences when the critical burning rate is exceeded, 
r^>r^^, such that subsurface melting can no longer propagate. At this time step, 
is set to zero and the existing melt layer thickness is allowed to decompose 
away in accordance with the time integral of r^ in succeeding time steps. This 
stage is completed when the melt layer has disappeared. 

The formation of the craters impacts both r and in accordance with the 
steady-state model. However, the development is allov/ed to proceed continuously 
with time (in the steady-state model it appears as a discontinuity with pressure). 
The changing surface area in the course of crater development is expressed as 
follows: 


''b 


AbF ^ "frozen", or cratered surface burn area 

A|j = geometric, or input burn area 

(S /S ) = normalized oxidizer surface exposed for a planar 

M surface, from Ref. (7). 

(S /S ) F = normalized oxidizer surface exposed for a 
cratered surface, from Ref. (7), 


(9) 
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The parameter (Sqj^/Sq)[: changes with time because the oxidizer (nitramine) 
penetration depth, h, changes with time. The relationship between (Sqj^/Sq)p 
and h appears in Ref. (1). The variation of the penetration depth with time 
may be expressed as follows: 

dh = (r^j^ - r^) dt (10) 

•■ox = V (-^ox/V,) (”> 

r-f = Af exp (-Ef/R„Tj) (12) 

r = nitramine monopropellant burning rate 

r^ = binder regression rate 

A = kinetics prefactor for the nitramine 

U A 

E = activation energy for the nitramine 

Ua 

= kinetics prefactor for binder 
= activation energy for binder 

A limit is placed upon the extent of penetration in accordance with the steady- 
state model. This limit results in a maximum or 

the steady-state model, it is known that the area change from to 

(S /S )p is the essential manifestation of the break point wherein r jumps 
from a value influenced by the binder to a value close to the nitramine mono- 
propellant rate. Therefore, in the course of the break point transition, it 
may be assumed that F is proportional to the change in (S.„/S ). This assumption 
is made in order to avoid jump discontinuities in the time-dependent model, and 
allow smooth (continuous) transition. Eq. (2), as applied during this third 
stage of surface development, becomes: 


ri 

"m 






r— 1 

[(-Ee/R) (1/T^ - 1/Tjl - 


1 


^^ox^F ^^ox^M 


L ^ ^ ^ i 




i^ox^F"^^6x^M_ 

) 


(13) 


rp= post-break steady-state burn rate 
rj^= pre-break steady-state burn rate 
(^^)p ~ value associated with Fp 
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It is anticipated that the break point in the time-dependent model need not jl; 

occur at the same pressure as in steady-state. That is, the-pressures for 

jf :■ 

r.>r and F> r„ may be different. Thus it is necessary to include values of p 

1 m m 

F which extrapolate planar surface values (F^;|) to higher pressure and cratered > 

surface values (Fp) to lower pressure. If the cratering actually occurs at .v 

some different pressure than steady-state, a proper value of F may still be used 
in accordance with the surface existing at the time. 

L'i: 

The additional unknowns presented by nitramine propellants are F 

r^, h, (Sqj^)p and A^^p. Eqs. (9-12) and two equations from the steady-state 
model are included in the forward difference scheme to provide the relationships. ;• 

Eq. (13) substitutes for Eq. (2), and Aj^p is used in Eq. (6). The changing surface [1 

structure is an added means by which the instantaneous burn rate may differ from f 

the equilibrium value, and a difference between Aj^ and Aj^p would bear upon effective 
burning rates deduced from the closed vessel experinent. The effective burning | 

rate is obtained by multiplication of Eqs. (9) and (13). ■: 


3.3,4 Computer Programming 

A computer program providing for the numerical solution of the equations was 
v/ritten for operation on a UNIVAC 1108 computer. The program has been successfully 
rendered operational and applied to test cases. An important limitation is the 
run time which can be encountered at very high burning rates, because a rapid 
transient demands a smal 1 time step for accuracy. To minimize the impact of this 
dependence, a scheme which automatically adjusts the time step and the solid 
propellant grid to the existing transient was successfully developed. The run 
time for a case in which the burning rate reached 20 in/sec was 15 minutes. A 
description of the computer program is presented in Appendix B. A source deck 
and listing can be furnished upon request. 
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•3.4 MODEL RESULTS 


3.4.1 Homogeneous Propellant Test Case 

The homogeneous propellant test case selected for model validation is "Test 
Propellant A" as specified by the JANNAF Workshop on Burning Rate Measurements, 
Memoranda of 26 January 1977 and 5 April 1977. This case was selected because of 
the extensive data provided in association with the purposes of the Workshop. 

Burning rates v/ere reported by six different laboratories for a given set of 
specifications and closed vessel pressure-time data. However, the identity of 
the propellant was not disclosed and strand burning rates were not given. In the 
absence of strand data, the reported burning rates were used as baseline equilibrium 
data for purposes of model application. The data were reported over the pressure 
range 5Kpsi-40Kpsi . An average of the results for the six laboratories was used at 
each pressure, and an extrapolation was performed to complete the input down to 
ambient pressure. Although these data are closed vessel data, the model would 
still be useful to show the effect of the closed vessel environment, if any. If 
such effect is small, as expected, the model should reproduce these data as well 
as the measured pressure-time curve. 

The burning rate results are shown in Figure 3. The solid line represents the 
actual closed vessel data, except for pressures below 5Kpsi wherein it is an 
extrapolation. The model prediction is represented by the circles. It is observed 
that the prediction is very close to the line between BKpsi and 40Kpsi, indicating 
a very small effect of the transient environment. In fact, the effect of the 
transient environment disappears completely for this test case above 13Kpsi. 

An interesting undulation appears in the predicted results between 4Kpsi and 
12Kpsi. It is not certain whether this reflects a real mechanism, a tradeoff 
between the transient combustion zone he'.ting and transpiration cooling during 
the pressure rise, or is simply a result of the forward difference numerical method. 
Similar undulations have appeared in computerized data reduction for closed vessel 
burn rate (e.g.. Ref. 24). 
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Figure 3, Predicted Effect of Closed Vessel Environment on Burn Rate 
of a Homogeneous Propellant 





The most interesting behavior appears at low pressure. The very high initial 
burning rates are largely a consequence of the assumed igniter. The model does 
not include a representation of the igniter as such, but uses an arbitrary 
input heat flux to get the burning under way. If this heat flux remains operative 
during burning, and is large compared to the combustion-induced heat flux, the 
burning rate will be high. The igniter is turned off at 1000 psi, and the 
burning rate immediately drops because the transpiration cooling becomes excessive 
for the remaining heat imposed. The burning rate climbs again v/hen the transient 
heating is re-established. Here, then, is a second type of burning rate undulation. 

A third type {not shown here) results from pre-ignition thermal soak of the propellant, 
and is sometimes referred to as the ignition "over-shoot" (25). Thus, there are 
mechanistic bases for low pressure burning rate undulations, depending upon the 
nature of the ignition process and developing combustion-induced heating. These 
effects, although real, are undesirable for purposes of closed vessel data 
acquisition and are lumped into the area of uncertainties. As a result, burning 
rates deduced from the closed vessel generally exclude results below 5Kpsi. An 
accurate quantification of the low pressure behavior v;ould require a model of the 
igniter. 

Predicted pressure-time results are shown together with data in Figure 4. 

Time is measured from ignition time rather than zero time to remove the difference 
in ignition delay and afford a better visual comparison. Initial pressurization 
is overpredicted because of the effect of the assumed igniter. This also causes 
the predicted web to advance relative to the actual web at a given time and 
pressure. The curves then merge and the prediction becomes very good at the 
higher pressures. The pressure results are consistent with the burning rate results. 

3.4.2 Nitramine Propellant Test Case 

The nitramine propellant test case selected for model validation is the 
92-micron HMX propellant reported by Lenschitz (Ref. 20). Much of the required 
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input information is contained within that paper, including actual strand burning 
rate data to serve as the equilibrium baseline. The paper also reports the 
pressurization history for that test. Other information required was furnished 
by Lenschitz. 

The burning rate results are shown in Figure 5. Plotted are the strand data 
(solid line), deduced closed vessel data (solid circles), and the predicted actual 
(open circles) and effective (x) burning rates in the closed vessel. As expected, 
there is a significant difference between strand and closed vessel burning rates in 
the case of a nitramine propellant, and the model serves to explain why. 

First, let us examine the region below 2000 psi. Here, much the same be- 
havior appears as was discussed for the homogeneous propellant; the burning rate 
is first dominated by the behavior of the assumed igniter, and drops when the 
igniter is turned off. Note, iiowever, the appearance of a second burning rate 
undulation just below 2000 psi. This is a consequence of a fourth mechanism for 
undulation, and will be unique to nitramine propellants v^hich exhibit a break 
point. What is happening here is that a break point is being traversed, causing 
a greater burning rate than can momentarily be sustained by the combustion heating 
than in existence. The rate falls until sufficient heating develops to re-sustain 
it. If the rate fell to the extent that the surface melt layer would re-form, 
then an endotherm v/ould reappear to drop the rate further; such a process might 
appear as a fifth mechanism for undulation. No data are reported belov/ 2000 psi- 
but the model results are not expected to be accurate in that regime because the 
igniter is not modeled. 

The region between 2000 psi and 10 Kpsi is interesting because it shows how 
the predicted actual burning rates in the closed vessel differ from equilibrium 
data. Basically, the post-break-point surface structure cratering of the prooellant 
is developing continuously and gradually in the course of the transient in the 
closed vessel. This continuous development lags the step-wise cratering at the 
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discrete test-to-test pressure intervals in the strand burner. Thus, the severity 
of the break point, as measured by the exponent in this transition region, is 
lessened by the closed vessel environment. At pressures above lOKpsi, the pre- 
dicted actual burning rates in the closed vessel merge with the equilibrium strand 
data as the surface structure becomes fully developed. 

Comparison between predicted and measured burn rates in the closed vessel re- 
quires use of the predicted effective burn rate, because that parameter accounts 
for the effective propellant burn surface area due to the cratering. The predicted 
effective burn rate does not depart significant from the predicted actual burn rate 
until pressures in excess of 3Kpsi. Over the range 3Kpsi-20Kpsi , the predictions 
are generally on the low side, but agree rather well with the data. Thus, it is 
concluded that the high closed vessel burn rates at high pressure, relative to the 
strand data, are due entirely to a surface area effect and not a transient burn 
rate effect. Lenschitz had rendered a similar conclusion (Ref. 20), but thought 
that the propellant grains were breaking up. Grain breakup is not required to 
explain these results if the cratered surface structure is recognized. 

The predicted pressurization is compared with data in Figure 6. Here, the 
data v/ere reported as rate of change of pressure versus pressure rather than as 
pressure versus time. This is a more severe test of the model because a number of 
factors enter into the rate of pressure change at a given pressure. The prediction 
falls below the data, which is consistent with the underprediction of the effective 
burn rates, but the agreement is rather good considering all of the factors involved. 

Based upon the results of this effort and the steady-state modeling effort, 
it would be expected that closed vessel burning rates should agree with strand 
burning rates when the HMX particle size is sufficiently fine that there is no 
break point or surface cratering. 


RATE OF CHANGE OF PRESSURE, Mpsi/sec 



Figure 6, Comparison of Predicted and Experimental Pressurization for 
the Test Corresponding to Figure 5 





SECTION 4 


COMBUSTION INSTABILITY CHARACTERISTICS OF NITRAMINE PROPELLANTS 

4.1 BACKGROUND HYPOTHESIS 

A belief that HMX may have favorable effects upon combustion instability 

driving, relative to AP, is based upon the qualitative information presented in 

Table 1. With respect to some factors which may be important, HMX appears to 

be more similar to ammonium nitrate (AN) and potassium perchlorate (KP) than to 

AP. Propellants based upon AN or KP have not been reported to be unstable 

(Ref. 26). HMX is similar to AN and KP in that there is a substantial oxidizer 

1 

melt laye>^ and particle size effects upon burn ^ate are relatively weak. 

Hov/ever, HMX surpasses AP as an energetic monopropel\ant. To the extent that 
the products of combustion include particulates which may afford acoustic damping, 
HMX would tend to produce more carbonaceous residue than AN (but perhaps no. 
as effective as spherical KCl droplets). At the commencement of this work, 
there was no reported systematic study of the effect of HMX on combustion driving. 

4.2 PROPELLANT SELECTION 

Standard thermochemical calculations and burning rate model calculations 
were performed to provide selection of the propellant formulations for T-burner 
combustion instability experiments. The purpose is to determine the effect of 
replacing ammonium perchlorate (AP) v/ith HMX on the response function of inert 
binder (HTPB) propellants, and of HMX addition on the response function of 
nitrocellulose (NC) propellants. In making these comparisons, it v/as desired 
to maintain propellant energy and burning rate reasonably constant within each 
group. Although differences can be corrected in tl:e course of interpreting the 
response function, it is considered that excessive differences will render the 
comparison less meaningful due to mechanistic bias. 

1 True at pressures below \:he HMX break point. Significant changes occur above 
the break point, which v;ould alter the hypothesis. Pressures of interest in 
solid rockets, as opposed to guns, are generally below these break points. 
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T.V3LE 1 


PREMISES FOR COMBUSTION STx\BILin' OF HMK PROPELLAKIS 


Propensity Toviard Instability Yes 

Substantial Oxidizer Melt No 

Bum Rate as Function of Strong 

Particle Size 

Condensed Phase Heat Release Strong 

Monopropellant Flame Energy Weak 

Particulates in Combustion Zone No 


M 


ilMX 

No 

No 

9 

Yes 

Yes 

Yes 

Weak 

Weak 

Weak 

Wej.k 

None 

Strong 

Weak 

None 

Strong 

?* 

Yes 



* - Depends on fuel~richness . For a given solids loading, an HMX propellant 
would be more fuel-rich than an AN propellant. Although it is possible 
to formulate such that no free carbon would be predicted by equilibrium 
thermochemistry, residues are sometimes observed. 


2 

to 

nj O 

to^ 

^ > 

^ Q 

^ © 


AP and HMX produce different thermochemical effects in HTPB binder because 
AP is an oxidizer whereas HMX is an energetic monopropellant. HMX propellants 
tend to be low burning rate propellants, but significant burning rate control is 
available in AP propellants through particle size adjustment. Tradeoff study 
produced the tv/o groups of selected HTPB propellants listed in Table 2; 
a low burning rate group (XA) and a high buicning rate group (XB). It is difficult 
to match AP propellant burning rates with HMX alone, so the burning rate of 
formulation XA-3 is lower than desired for its group although its energy level 
is adequate. The two mixed oxidizer propellants are more closely related to their 
AP analogs in burning rate and energy. 

It was found that the addition of meaningful quantities of HMX to an inert- 
plasticized MC binder produced considerable increases in propellant energy. Use 
of a relatively stable energetic plasticizer (TMETN) reduced the impact of 
HMX addition on propellant energy, and was considered more acceptable for 
research purposes than nitroglycerine. The resulting propellants comprise the 
third group (XC) listed in Table 2. These propellants v/ere not tested in the 
course of this program because additional testing of the HTPB propellants was 
found to be necessary in order to acquire sufficient meaningful data. 

4.3 EXPERIMENTAL RESULTS 

Stability tests were carried out in 2 1/2-inch diameter (I.D.) T-burners 
coupled directly to a surge tank. The propellant configuration used was a 3/8- 
inch thick disk. Data v^ere obtained using the growth-decay technique (Ref. 27). 

For each propellant, tests were conducted at nominal frequencies of 500, 900, 
and 1900 Hz anu at pressures of 500 and 1000 Psi. Several tests were performed 
at each condition. Many of the tests did not produce oscillations, or the data 
did not conform with T-burner standards, and those tests are excluded from the results. 


26 


TABLE 2 


wt-: 


PROPELLANT FORMULATIONS FOR T-BURNER TESTS 



XA-1 

XA-2 

XA-3 

XB-4 

XB-5 

XC-6 

XC-7 

AP(9iO 


- 

- 

41.25 

31.25 

- 


AP(50y) 

40.00 

20.00 


- 


- 

“ 

AP(90p) 

- 

- 

- 

41.25 

31.25 

- 

- 

AP(200p) 

40.00 

20.00 

- 


- 

- 


HMX(Ap) 

- 

20.00 

40.00 

- 

22.50 

- 

50.00 

HMX(15p) 

- 

20.00 

40,00 

- 

- 

_ 

- 

R-A5/HTPB 

20.00 

20.00 

20.00 

17.50 

15.00 

- 

- 

NC/TMETN 

" 

- 

- 

- 

- 

100.00 

50.00 

(T„/M)^^^ 

10.34 

10.30 

10.32 

10.70 

10.77 

10.49 

10.86 


0.26^ 

0.21^ 

0.16^ 

0.48^ 

0.43^ 

0.30 

0.33 


1 — Measured Values 


The results for the lov/ rate and hi^n rate HTPB propellants are given in 
Figures 7 and 8, respectively. The results are plotted as the sum of the measured 
growth and decay coefficients vs. frequency. The propellants v/ere all relatively 
stable in that the measured growth coefficients exceeded the decay coefficients 
only for the two high rate propellants at the 1900 Hz, 500 psi test condition. 

Such behavior may subject the absolute results to criticism, but the relative 
results are considered valid. 

The low rate XA-1 propellant was the weaker acoustic driver of the two 
propellants which did not contain HMX. The low rate results from the use of 
coarser particle sizes. Its results appear to be relatively independent of 
pressure. The minimum between 500 Hz and 1800 Hz is considered real because only 
one test out of seven conducted at 900 Hz produced measurable acoustic driving, 
whereas most of the excluded data at the other frequencies stemmed from poor 
raw data. 

The two low rate propellants containing HMX exhibited unusual results. At 
500 Hz, the pressure traces showed an immediate grov/th of oscillations following 
ignition to amplitudes of about 2 % of the mean and then a slow decay over the 
duration of the burn. At 1900 Hz, the tests exhibited no oscillations at all. 

At 900 Hz, the 500 psi tests were like those at 500 Hz and the 1000 psi tests 
were like those at 1900 Hz. 

Post-test examination of the XA-3 tests showed the T-burner to be full of a 
layered, carbonaceous char. In the 1900 Hz tests, this char appeardd to have 
been pulverized. For the lower HMX concentration XA-2 propellant, the char v;as 
less extensive, merely coating the wall of the T-burners rather than filling the 
interior. The char is presumed to be due to the fuel -rich nature of the two low 
rate propellants containing HMX. It was somewhat surprising because prior 
experience with similar propellants did not encounter char. A review of the 
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formulations revealed that the prior propellants were much more heavily plasticized 
gun propellants, whereas the propellants formulated here were more consistent 
with the structural integrity requirements of rockets. Another consequence of 
the reduced plasticizer, or the char, was that the burning rate was about 20% 
lower than expected. The char is indicative of non-equilibrium thermochemistry 
in that there should be no free carbon at a solids loading of 80% HMX. In any 
event, the contribution of HMX to stability could not be identified between 
combustion driving or mechanical suppression from these tests. 

The high rate propellant data in Fig. 8 is mo^’e conclusive as to the role 
of HMX. These propellants did not produce char. Both XB-4 and XB-5 showed 
clear pressure and frequency effects, increasing in driving strength with de- 
creasing pressure and increasing frequency. The XB-5 HMX-containing propellant 
would not drive oscillations at all at the 1000 psi - 500 Hz test condition, 
and its data at 500 psi and 1000 psi fall below the respective data for the 
all-AP oxidizer XB-4 propellant. The measured damping coefficients at each 
frequency v/ere about equal for the two propellants, indicating that the differences 
in acoustic driving are real. More data are required for clean-burning propel- 
lants, and testing of the XC group is included in future plans. 


30 


SECTION 5 


EFFECT OF NITRAMINE INGREDIENTS ON STEADY-STATE BURNING RATE 
(Data Supplement to Annual Progress Report AF0SR-TR-76-1136) 

5.1 PREFACE 

Strand burning rate data supplementing the data reported In the FY 1976 
Annual Progress Report (Ref. 1) were furnished by AFATL toI lowing the publica- 
tion of that Report. These data are presented and discussed in the following 
subsections. It is concluded that the data provide further confirmation of the 
results and predictions cf the steady-state model developed in the course of 
that research. 

5.2 EFFECT OF HMX PARTICLE SIZE IN ACTIVE BINDER PROPELLANTS 

Figure 5 of Ref. (1) showed the effect of HMX particle size on the burning 
rate of nitroplasticized polyurethane (NPPU) propellant. Analogous data for NC 
propellant are now available and are presented in Fig. 9. The same basic trends 
appear as in the NPPU propellant, as was expected. The fine propellant exhibits 
a relatively uniform pressure exponent over two orders of magnitude of pressure. 

The coarse propellant begins to depart significantly at a pressure above IKpsi, 
with an upward break in pressure exponent covering the range 1.3Kpsi to 2.5Kpsi. 

A second upward break appears to begin at 4Kpsi and continue to high pressure. 

As a result, the burning rates of the coarse propellant have a multiple ziq-zaq 
appearance and are higher at high pressure. The burning rates of the NC propellants 
are higher than those of their NPPU analogs, vyhich is attributed to the higher 
energy and burning rate of the NC binder. 

A similar trend is shown by the mixed HMX-TAGM propellants in Fig. 10. A 
portion of the data for the propellant containing the fine HMX appeared in Fig. 7 
of Ref. (1). No data were available for the propellant containing the coarse 
HMX. Results for both propellants are now complete. The propellant 
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PRESSURE, PSl 


Figure 10 Effect of HMX Particle Size in Mixed KMX-TAGN, NC Propellants 






containing fine HMX exhibits a relatively uniform shape over a wide pressure 
range. The propellant containing coarse HMX exhibits an upward break in pres- 
sure exponent covering the range IKpsi to 3Kpsi, such that its burning rates are 
higher at high pressures, as was expected. Comparing Figs. 9 and 10, the effect 
of the TAGN is to raise burning rates, more so at low pressure than at high 
pressure, such that the overall net pressure exponent is reduced. This effect 
has been attributed to the faster kinetics of TAGN decomposition and combustion 
as a monopropellant (Ref. 1). 

Figs. 9 and 10 illustrate the effectiveness of fine particle size as an 
approach to achieve desirable slope uniformity. Fig. 9 confirms that it is 
possible to do so with HMX alone if the particle size is sufficiently fine. 

Fig. 10 confirms that TAGN addition is helpful provided that the HMX component 
does not exhibit a slope break. 

5.3 EFFECT OF TAGN IN ACTIVE BINDER PROPELLANTS 

Data presented in Figs. 9 and 10 allow the completion of Fig. 7 of Ref (1). 
This completed figure is presented as Fig. 11. The effect of TAGN substitution 
for HMX, in increasing amounts, is thereby illustrated. It v/as predicted that 
the data for the 29% TAGN propellant would interweave v/ith the data for the 
50% TAGN propellant at the higher pressures. This prediction is verified by 
the data. A break point appears for the 50% TAGN propellant at lOKpsi. This 
was explained by the burning rates reaching such a level as to trigger the break 
in the HMX component, and was described as a disadvantage of excessive TAGN 
addition. In view of that explanation, the 29% TAGN propella.it also should 
exhibit a break point when the burning rates reach the same level. This is not 
so obvious in its high pressure data, perhaps because the region is not as well 
characterized with data points. In any event, it is confirmed that TAGN addition 
raises burning rates, and that optimum amounts are helpful as an approach to 
achieve desired slope uniformity. 
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There has been a question as to whether or not propellants containing 
alone v/ould exhibit a break point. The uppermost data in Fig. 11 are for such 
a TA6N propellant. The model predicted that a break point exists near the region 
where considerable data scatter appears, and the scatter was interpreted to be a 
manifestation of the break. The model also predicted that the second upward 
break in the 50% TAGN propellant is due to the TAGN. The question, however, 
could not be resolved from the data. Complete data are now available for the 
'■ 3 % TAGN, NPPU propellant (cf. Fig. 6 of Ref [1]). These data are shov/n to- 
gether with the Fig. 11 data in Fig. 12. A break point for this TAGN propellant 
clearly appears at lOKpsi. Therefore, it can be concluded that TAGN propellants 
can exhibit break points. Ironically, this break point can be mitigated by 
lowering the burning rates, through HMX substitution, which would yield the GAU-S 
propellant. The shape of the NPPU propellant burning rate curve lends credence 
to the interpretation of a high slope (break) region in the case of the NC 
propellant. The NC propellant has higher burning rates because of the greater 
TAGN concentration and the more energetic binder. 


5.4 MATCHING NITRAMINE AND ACTIVE BINDER 

It was pointed out in Ref. (1) that matching the burning rates of the active 
binder and nitramine as monopropellants could yield a propellant having a uniform 
pressure exponent. The rationale was that a break point would have no way to 
manifest itself, and so would not appear in the burning rate data. In support 
of this hypothesis, data and results of model calculations were presented in 
Figs. 31-33 of that Report. These figures showed trends as the nitramine 
burning rate approaches and falls below the binder burning rate. It was suggested 
that a propellant combining fine EDNA and a tailored active binder having the 
same burning rates as EDNA would be a practical test of the approach. 
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5^ TAGN 




O - 75% TAGN/NC 
@ - 63% TAGN/NPPU 



PRESSURE, PS) 


Figure 12 • Burning Rates of TAGN Active Binder Propellants 
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Unfortunately, the propellant that was furnished combined coarse EDNA and 
the baseline NC binder. As a result, the binder burning rates are lower than 
the EDNA burning rates in accordance with the Figure 31, and the break point 
mechanism would be triggered at a low pressure. With this propellant, however, 
it is possible to complete the EDNA propellant portion of Fig. 31 ard these data 
are shown in Fig. 13. The premise that the post-break propellant burning rates 
follow the nitramine monopropellant burning rates is confirmed. Without data at 
lower pressures, v/here. a" break would be exoected, it does not aooear to matter 

t . 

that the binder but'ning rates are low. A transition or jump upward to the EDNA 
line is not seen over this pressure range, such that this EDNA/NC propellant 
appears to achieve a desirable result. 

— Aft illustration of binder tailorability to implement the foregoing approach 
is presented in Fig. 14. The lov^er dashed line is for nitrocellulose (NC) 
containing an inert plasticizer. The upper dashed line is for nitrocellulose 
containing a very energetic plasticizer. The line fitting AFATL data is for 
nitrocellulose containing an energetic (intermediate) plasticizer. The range of 
burning rates is observed to cover a factor of 3. A larger range would be available 
by adjusting ingredient proportions, subject to energy and processing limitations; 
a practical range of 5 is indicated by existing active binders. The intermediate 
binder shown in Fig. 14 has burning rates quite close to EDNA burning rates, and 
v/ould be an attractive candidate. Note, however, that the details of the binder 
data indicate some break point mechanism in the binder which might be reflected 
by the propellant using it. It is difficult to comment upon break points in active 
bindiTS in general because the data are usually reported as straight lines. 
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Figure 13, 


Comparison of Monopropellant Burning Rates and EDNA 
Propellant Burning Rates 
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SECTION 6 


PRESENTATIONS AND PUBLICATIONS 


The following presentations and publications have been generated by this 
research contract: 

(1) Cohen, N.S. and Strand, L.D., "Nitramine Propellant 
Research", 13th JANNAF Combustion Meeting (CPIA Pub- 
lication ?81 , Vol. I, Dec. 19'76) pp. 75-87. 

(2) Cohen, N.S. and St»'and, L.D., "Nitramine Smokeless 
Propellant Research", 1977 AFOSR/AFRPL Rocket 
Propulsion Research Meeting, L^^ncaster, CA (Mar, 1977). 

(3) Cohen, N.S., Strand, L.D., and Price, C.F., "Analytical 
Model of the Combustion of Multicomponent Solid Propellants", 
AIAA Paper 77-927, AIAA/SAE 13th Propulsion Conference, 
Orlando, FL (July, 1977). 

(4) Cohen, N.S. and Strand, L.D., "Transient Combustion 
Response Characteristics of Nitramine Smokeless 
Propellants", 14th JANNAF Combustion Meeting , Colorado 
Springs, CO (Aug., 1977), publication pending. 


SECTION 7 


CONCLUSIONS 


A comparison of strand and closed vessel burning rate data for homogeneous 
propellants reveals that effects of pressurization are minor, and that the two 
methods should give consistent results at pressures following the ignition 
transient in the closed vessel. A description of burning rates during the 
ignition transient in a closed vessel would require a model of the igniter. 

The burning rates of nitramine propellants as measured in a closed vessel 
will differ from strand data for those propellants which exhibit a break point. 

The effect of the transient environment in a closed vessel is to create a time 
lag in the development of the post-break-point cratered surface structure, which 
results in a less severe exponent shift than shown by the strand data. The 
reduced severity of the exponent shift will have a practical benefit in transient 
devices. However, the cratered surface structure will maintain high mass burning 
rates and exponents because of the increased burn surface area. Methods to avoid 
or mitigate the cratering, as discussed in the Annual Progress Report for FY 1976 
{Ref. 1), are still recommended. Additional data have been presented here which 
confirm the conclusions of that Report. 

Several mechanisms have been identified v/hich might explain burning rate 
undulations often observed from closed vessel data reduction, in particular at 
low pressure: ignition overshoot, termination of the igniter, tradeoff between 

developing combustion zone-induced heating and increasing surface transpiration 
plus, in the case of nitramine propellants, break-point overshoot and melt layer 
refo'^mation. 

HMX-containing composite propellants are more stable than comparable AP 

propellants. One mechanism may be related to mechanical suppression due to the 
* 

generation of solid material combustion products in certain fuel-rich formulations. 
However, it also is true that HMX reduces the combustion driving in clean-burning 
formulations tested thus far. Additional study is recommended. 


SECTION 8 


REFERENCES 


1. Cohen, N.S. and Strand, L.D., "Nitramine Propellant Research", Report AFOSR- 
TR-76-1163, JPL TM 33-801, Jet Propulsion Laboratory, Pasadena, CA (Oct. 
1976); also cf. AIAA Paper 77-927 (July, 1977). 

2. Heiney, O.K. , "Advanced Gun Propellants", J. Nat. Defense, Weapons Tech- 
nology, pp. 152-157 (Sept. - Oct. 1973). 

3. "Research Objectives Document", HQ AFSC TR 76-01 (May 1976). 

4. Meier, T.C., "Air Force Rocket Combustion Research", Presentation at 13th 
JANNAF Combustion Meeting, Monterey, CA (Sept. 1976). 

5. Steifel, L., "Reviev^ of Workshop on the Combustion of Nitramine Propellants 
for Guns", 10th JANNAF Combustion Meeting (CPIA Publication 243, Vol. I, 

Dec. 1973), pp. 199-216. 

6. Kumar, R.N. and Strand, L.D., "Theoretical Combustion Modeling Study of 
Nitramine Propellants", J. Spacecraft and Rockets 1^, 427-433 (July, 1977). 

7. Cohen, N.S. and Price, C.F., "Combustion of Nitramine Propellants", J. 
Spacecraft and Rockets J2_, 608-612 (Oct. 1975). 

8. Simmons, R.L., "Workshop Report: Nitramine Gun Propellant Combustion", 

13th JANNAF Combustion Meeting (CPIA Publication No. 281, Vol. I, Dec. 

1976) pp. 1-8. 

9. "HMX Propellant Combustion Studies", Contract F04611-75-C-0019, Hercules, 
Inc., Magna, UT (in progress). 

10. "Improved Gun Propellant", Contract F08635-76-C-0099, Rocketdyne Division 
of Rockwell /International Corp., Canoga Park, CA (in progress). 

11. Heiney, O.K. , Presentation at the AFOSR-AFATL Explosive Combustion Meeting, 
Eglin AFB, FL (Apr. 1974). 

12. Cohen, N.S., "Report of Workshop on Combustion Instability of Smokeless 
Propellants", 14th JANNAF Combustion Meeting (to be published). 

13. Von Elbe, 6., "Theory of Solid Propellant Ignition and Response to Pressure 
Transients", 19th ICRPG Solid Propulsion Conference (CPIA Publication 18, 
Vol. Ill, 1973), pp. 157-181. 

14. Coates, R.L. and Cohen, N.S., "Correlation of Solid Propellant Extinguish- 
nient Data", 3rd ICRPG Combustion Conference (CPIA Publication 138, 

Vol. I, Feb. 1967), pp. 213-218. 

15. Krier, H., T'ien, J., Sirignano, W. and Summerfield, M., "Nonsteady Burning 
Phenomena of Solid Propellants", J. AIAA 278-285 (Feb. 1968). 


43 


16. Cohen, N.S., "Ballistic Predictions for Mass -Augmented Solid Rocket Motors", 
Report AFRPL-TR-71-133, Lockheed Propulsion Company (Dec. 1971). 

17. Cohen, N.S, and Harry, D.P. Ill, "Hydroquench Thrust Termination Analysis", 
Report AFRPL-TR-73-62, Lockheed Propulsion Company (Aug. 1973). 

18. JANNAF Workshop on Burning Rate Measurements, Calculations and Data Reduction 
Procedures, Ballistics Research Laboratories, Aberdeen Proving Ground, MD 
(Oct. 1976). 

19. Rocchio, J.J., et al., "Results of Recent Experimental and Theoretical Studies 
of Nitramine Gun Propellant Performance", ibid. Ref. (8), pp. 9-19. 

20. Lenschitz, C., "The Role of Thermochemistry in HMX Propellant Burning", 

12th JANNAF Combustion Meeting (CPIA Publication 273, Vol . II, Dec. 1975), 
pp. 301-321. 

21. Derr, R.L., Presentation at the AFOSR-AFATL Explosive Combustion Meeting, 

Eglin AFB, FL (Apr. 1974), 

22. Price, C.F., private communications, Ballistics Research Laboratories, 

Aberdeen Proving Ground, MD (Mar. 1976). 

23. Grollman, B.B. and Nelson, C.W., "Burning Rates of Standard Army Propellants 
in Strand Burner and Closed Vessel Tests", ibid. Ref. (8), pp. 21-43. 

24. Mitchell, S.E. and Horst, A. VI. Jr., "Comparative Burning Rate Study", ibid 
Ref. (8), pp. 383-404. 

25. Kooker, D.E. and Nelson, C.W., "Numerical Solution of Three Solid Propellant 
Combustion Models During a Gun Pressure Transient", 12th JANNAF Combustion 
Meeting (CPIA Publication 273, Vol, I, D'=«c. 1975), pp. 173-197. 

26. Price, E.W. and Culick, F.E.C., "Combustion of Solid Rocket Propellants", 

AIAA Professional Study Series, p, 242. 

27. "T-Burner Manual", CPIA Publication 191 (Nov. 1969). 





44 


APPENDIX A 


NOMENCLATURE 


A 

A, 


bF 


c 


"ox 

^ext 


h 

h 

c 

k 

n 

P 

* 

P 

V 

Qc 


;ss 


r 

r. 


"i 


ox 


R 

^^ox^ F 
t 

T 


kinetics prefactor for propellant decomposition 
propellant burn surface area 

nitramine propellant burn surface area following cratered surface development 

kinetics prefactor for binder decomposition 

kinetics prefactor for nitramine decomposition 

surface area of closed vessel wall exposed to heating 

propellant heat capacity 

gas heat capacity 

activation energy for propellant decomposition 
activation energy for binder decomposition 
activation energy for nitramine decomposition 
igniter heat flux 

depth of surface cratering in nitramine propellants 

heat transfer coefficient for v/all heatinq 

propellant thermal conductivity 

mass of gas in chamber 

burning rate pressure exponent 

pressure 

rate of change of pressure 

heat flux to propellant surface 

propellant heat of decomposition 

rate of heat loss to closed vessel wall 

steady-state burning rate at P 

regression rate of binder 

steady-state post break point burn rate at P 

instantaneous propellant burn rate 

steady-state pre-break point burn rate at P 

regression rate of nitramine 

gas constant 

universal gas constant 

fract' on of exposed nitramine surface on a planar melt surface 
fraction of exposed nitramine surface on a cratered surface 


value of associated with r^ 
time 

temperature 

A-1 
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propellant adiabatic flame temperature 

instantaneous gas temperature 

initial (ambient) temperature 

instantaneous propellant surface temperature 

steady-state propellant surface temperature associated v/ith 

closed vessel v/all temperature 

instantaneous gas volume in the clos''^ vessel 

propellant web 

distance from propellant surface into propellant 

radiation emissivity of gas 
covolume 

ratio of specific heats 
propellant thermal diffusivity 
propellant density 
Stefan-Boltzmann constant 


APPENDIX B 


COMPUTER PROGRAM DESCRIPTION 


B-1 SUBROUTINES 

The computer program file is identified as KT, and its map of linked 
elements for execution is identified as KRTT. The component subroutines 
are listed as follows: 


Segment 

MAIN PROGRAM 

Subroutine 

EVES 

Subroutine 

SETUP 

Subroutine 

DIFEQ 

Subroutine 

TSWALL 

Subroutine 

TSPROP 

Subroutine 

PRINT 

Subroutine 

FIND 


A schematic of the computer program elements is shox-m in Fig. B-1. 

Segment tIAIN is the master control element. This element reads input, 
computes certain parameters of an input nature, writes the input, and 
then transfer control to EVES . Subroutine EVES is a simultaneous 
differential equation solver that uses the Adams predictor -corrector 
method. EVES sets the time step, predicts ahead, compares the results 
and certifies accepted answers. Subroutine SETUP initializes all of the 
variables, and is the first subroutine called by EVES. Thus SETUP provides 
the initial conditions for the time-dependent problem. Subroutine DIFEO 
computes the time derivatives of gas pressure, gas volume, gas v/eight, 
gas temperature and distance burned. DIFEQ contains the time-dependent 
conservation equations which compute these differentials. EVES calls 
DIFEQ after the problem has been initialized. DIFEQ calls TSPROP to 
obtain the propellant combustion contributions to the conservation 
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a) MAX. RUN TIME EXCEEDED 
OR b) MAX. PRESSURE EXCEEDED 
OR c) WEB BURN-OUT ACHIEVED 


Figure B-1, Schematic of Computer Program Elements 
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equations, and TSWALL to obtain the heat loss to the wall which enters 
into Liiose equations. TSPRCP computes the time-dependent temperature 
profile in the solid propellant, the burning rate which is a function 
of surface temperature, and includes mechanisms by which nitramine behavior can 
influence the burning rate. A finite difference approximation of the 
transient heat equation is used to compute the temperature profile, which 
develops in response to igniter-induced heating and then to combustion 
heating. TSPROP may control the time step when high burn rate leads to 
steep temperatvire profiles. TSWALL computes the convective and radiation heating 
of the wall, and the temperature profile in the wall, in determining the 
heat loss. TSPROP and TSWALL return answers to DIFEQ, and DIFEQ returns 
answers to EVES. EVES then calls PRINT to update the accepted ansv;ers and to 
write them at pre-selected time intervals. Subroutine FIND is an auxiliary 
subroutine that is used to interpolate. 

The computations involving TSPROP, which is the essence of the 
combustion model, warrant a more detailed discussion which is provided in 
Section 2 of this Appendix. A source deck and listing of the computer 
program can be furnished upon request; it is liberally embellished with 
comment cards to assist the user. A manual for inputs is provided as 
Section 3 of this Appendix. Output is discussed in Section 4. 

B-2 THE COMBUSTION MODEL CALCULATIONS (TSPROP) 

The following sequential discussion describes the computations perfoinned 
in TSPROP. 

1. The time step is first determined by EVES. 

2. Tabular strand burning rates are initialized for the beginning of 
this time step. 
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3- The propellant web calculational mesh and counter are updated 
for the web which was burned away in the previous time step, 

4. A further adjustment to the mesh is made depending upon whether 
the reference mesh is greater than 1/10 the thermal wave thickness 
of the previous time step, or less than 1/25 of that thickness. 

Too coarse a mesh impairs accuracy, too fine a mesh will increase 
computer time by requiring a smaller time step. If the mesh is 
adjusted, the temperature profile is redefined by logarithmic 
interpolation to conform to the new mesh. This technique of 
adjusting the mesh has provided considerable savings in computer time, 

5. The time step required by the mesh is compared to the time step 
required by EVES, and the smaller value controls, 

6. The igniter "induced he fl’AX is gradually diminished to zero when 
sufficient burning has been achieved. The method of representing 
the igniter is arbitrary. There is no model of the igniter as 
such, only an input heat flux to start, 

7. The heat flux to the propellant is computed, consisting of igniter 
and combustic'n contributions. The combustion heat flux is equal to 
that V7hich would exist under equilibrium burning at the burning rate 
of the present time step, and therefore assumes a quasi-steady gas 
combustion zone. The appropriate tabular strand burning rate is 
updated to the present time step for this purpose. If the propellant 
is a nitramine propellant, and the actual burning rate of the 
previous time step is below the breakpoint burning rate, then the 
exothermic surface heat release term is modified by the heat of 
fusion of the nitramine. 
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Using a finite-difference expression j. the resulting heat flux 
is converted into an artificial surface temperature that is 
used to drive the thermal profile in the depth of the solid* 

9. The finite difference form of the transient heat equation is used 
to compute the temperature at each mesh location in the depth of 
the solid; including the true surface temperature. The rate of 
change of temperature is computed and multiplied by the time step, 

10. The actual burning rate is computed from the surface temperature 
using an Anhenius expression. 

11* Propellant density is known from inputs. The present burn area 
is kno™ from subroutine DIFEO. Knowing the actual burn rate 
the propellant weight flow rate is then computed, 

12. If the propellant is homogeneous, this burn rate and weight 

flow rate are returned to DIFEO at this point. If the propellant 
is a nitramine propellant, additional computations are performed. 

13. A test is made comparing the actual burning rate vjith the break -point 
burning rate. If the actual burning rate is less, the melt layer 
thickness of the nitramine is determined by the location of the 

melt temperature in the thermal profile* If the actual burning 
rate is now less but was at one time greater than the break point 
rate 5 such that the surface is now cratered, then an additional 
computation is performed to begin or continue to knit the surface 
back to a planar melt* This latter condition can result from de- 
pressurization or reductions in heat flux* 

14. After computing the melt layer thickness, the program returns to 
DIFEQ if the surface is now or Is once again, a planar melt* If the 
surface has not yet knitted completely, then adjustments are made 
according to the extent of cratering as discussed below. 


B-5 


15. If the actual burning rate is greater than the break-point 
burning rate in step 13, then the melt layer thickness is 
reduced by the product of the burning rate and the time step. 

If the melt layer thickness is still greater than zero, the program 
returns to DIFEQ. If the melt layer has disappeared by this 
mechanism, then the extent of cratering is computed according 
to the difference between the nitramine and binder burning rates. 
These ingredient burning rates are calculated from Arrhenius 
expressions for the present surface temperature. If the binder 
is an active binder, its rate is determined from the updated 
tabular input. 

16. If the surface is cratered, either as a result of step 15 or 
step 14, the actual burning rate is upgraded in accordance with 
the mechanism coFitained in the steady-state model. The 
parameters which produce this upgrade are computed from the 
known extent of cratering and t^ bular burn rates as presently 
updated. Associated v;ith this upgraded burn rate is an upgraded 
burn area, also computed from the cratering. Thus an upgraded 
weight flow rate is computed, and these upgraded values are 
returned to DIFEQ instead of the values from Step 12. 

B-3 CARD INPUT SEQUENCE AND STANDARD VALUES 

The following describes input cards in sequence, with appropriate 
comments and numerical values v/hich have been used. 


OMGINAEmEa 

OF BOOB QOAIJK 
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JOB CARD 

SEQUENCE OF CONTROL CARDS 
EXECUTION CARD 


CARD 


COLUMN 


ENTRY 


2 

3 


3 

2 

12 

22 

32 


42 


52 


32 

42 

52 

62 

2 

12 

22 

32 


SYMBOL 

ICASE 


UNITS 


Type of Propellant 

Comment: ICASE - 1 for homogeneous propellants 

ICASE = 2 for nitramine propellants 


Title of Case or Group 
Initial Chamber Gas Volume 
Initial Pressure 
Initial Temperature 
Igniter Heat Flux 


ITIT 

VINIT 

PINIT 

TINIT 

FEXT 


Comment: Values of 100-200 have been used. 

Reference Calculational Grid 
size of Propellant DELXO 

Comment: A value of 0.0002 

has been used, 

Calculational Grid Size for Wall DELR 

Comment: Divide v/all thickness 

by 10 


3 

62 

Chamber Diameter 

DIA 

4 

2 

Surface Area of Chamber Wall 

AWALL 


12 

Thermal Diffusivity of Wall 

ALPHAN 


22 

Flame Temperature of Propellant 




at 1000 psi 

TFP 


Comment: Divide isocharic flame 
temperature by GAMMA 
Gas Molecular Weight 
Ratio of specific heats 
Specific heat of gas at 10,000 
psi 

Specific heat of solid propellant 


Mt-T 

GAMMA 

CGP 

CSP 

COVOL 


2 

Ib/in 

2 

cal/cm -sec 


3_n, 


in* 


in.. 


ig. 

cm /sec. 


gm/mol 


cal 

o— 

gm- K 

cal 


-°K 


gm 
in /lb 


Microns 


Covolume 

Weight Fraction of Nitramine ALFA 

Mean Particle Size of Nitramine DAVE 

Maximum Surface Penetration 
Associated with Post-Break-Point 
Surface Cratering HDMAX 

Comment; Use a value of 3. for sizes greater than 100 microns 
Use cl value of 1. for sizes less than 10 microns 
Use a value of 2, for intermediate sizes ^ 

gm/cm^ 
gm/cm 
inVsec. 


42 

Nitramine density 

RHOX 

52 

Binder density 

RtiOF 

62 

Propellant Thermal Diffusivity 

ALPHAP 
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CARD 


COLUMN 


ENTRY 


SYMBOL 


UNITS 


6 


2 


12 


22 


32 


42 


52 


62 


7 


2 






12 




Ot 


22 


32 


42 


Reference Propellant Surface Temperature 

Temperature TSREF K 

Comment: For homogeneous propellants ; select 

a value consistent with the burning rate 
at 1000 psi and the kinetics expression 
for surface decomposition; an example is 
523‘^K for 0.35 in/sec. and an activation 
energy of 50000 cal/mol. Surface temperatures 
for nomogeneous propellants are in this 
vicinity. For nitramine propellants^ do 
the same but with a pre -break-point burning 
rate; for example, 926^K for 0,13 in/sec and 
50000 cal/mol. Nitramine surface temperatures 
are circa y00°K, and may be acquired from the 
steady-state model. 

Reference burning rate associated with RREF in/sec. 

TSREF 


Comment: See comment for TSREF 



Nitramine Melting Point 

Comment: The value for HMX is 551. 

The value for RDX is 467. 

TMELT 

°K 

Activation Energy for Nitramine or Homo 

- 


geneous Propellant Decomposition 
Comment: A value of 50000. has been 

used for each. 

Activation Energy for Binder 

ESP 

cal/mol 

De compo s it ion 

Comment: Use a value of 16900. for 

HTPB. For active binders, 
input 0. 

EF 

cal/mol 

Prefactor for Nitramine Decomposition 

Comment: A value of 0.5000 E f 10 has 
been used for HMX and RDX. 

ASP 

cm"^-sec 

Prefactor for binder decomposition 

AF 

gm 

2 

cm -sec. 

Comment: A value of 299 has been used 

HTPB. For active binders, input 0. 

for 

cal/gm. 

Surface Heat Release 
Comment: A value of -225. has been 
used to represent HMX and energetic 
homogeneous propellants. 

QSURP 

Nitramine Heat of Fusion 
Comment: A value of 132.5 is used 

OLM 

cal/gm 

for HMX and RDX, 


2 

Nitramine Thermal Diffusivity 
Comments: A value of 0,0011 has 

XKAPO 

cm /sec 

been used. 


2 , 

Binder Thermal Diffusivity 
Comment: A value of 0.0011 has 

XKAPF 

cm /sec 

been used. 

Nitramine Thermal Conductivity 
Comment: A value of 0.00049 has been 

used. 

XLAMO 

cal 

cm-sec- 
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CARD 


7 Cont*d 


9 


NEXT 


NEXT 


NEXT 


NEXT 


NEXT 


COLUMN 


ENTRY 


SYMBOL 


52 Binder Thermal Conductivity XLAMF 

Comment: A value of 0.00044 

has been used. 

62 Wall Thermal Conductivity XKN 


2 

2 

12 


Number of Points in the Table of 

Burn Area vs Web NWEB 

Web Burned 

Comment: The first table entry WEBlX 

should be zero 

Burn area ABX 

Comment: This is the total burn 
area at zero web burned. 


REPEAT FOR mEh CARDS TO DEFINE AREA VS WEB. 

2 Number of Points in the Table of NPl 

Pre-Break-Point Burning Rate 
vs. Pressure 

2 Pressure PMX 

12 Pre-Break Point Strand Burn Rate RMX 

Comment: Cover the full pressure 
range with the series of cards, 
extrapolating to PINIT and to the 
maximum pressure expected. The 
series may be omitted for homogeneous 
propellants, whence the NPl card 
should be omitted also. 


UNITS 

cal 

cm-sec- K 
cal 

o 

cm-sec“ K 


in. 


. 2 
rn 


Ib/in. ^ 
in/sec. 


REPEAT FOR NPl CARDS TO DEFINE PRE -BREAK-POINT 
BURNING RATE VS PRESSURE 

2 Number of Points in the Table of 

Post-Break-Point Burning Rate vs 



Pressure 

NP2 

— 

2 

Pressure 

PFX 

Ib/in. 

12 

Post-Break -Point Strand Burn Rate 

RPX 

in/sec 


Comment: Cover the full pressure 
range v;ith this series of cards, 
extrapolating to PINIT and to the 
maximum pressure expected. This 
series may be omitted for homogeneous 
propellants, v/hence the N?2 card 
should be omitted also. 

ORIGINAE 

PAGE Ib 

REPEAT 

FOR NP2 CARDS TO DEFINE POST-BREAK-POINT 

OP POOR QUALITY 

BURNING 

RATE VS PRESSURE 




2 Number of Points in the Table of 

Active Binder Burning Rate vs 
Pressure NP3 




CARD 

COLUMN 

ENTRY 

SYMBOL 

UNITS 


2 

Pressure 

PBX 

Ib/in^ 


12 

Active Binder Strand Burn Rate 

RBX 

in/sec 


Comment: Cover the full pressure range with this 
series of cards, extrapolating to PINIT and to 
the maximum pressure expected. This series may 
be omitted for inert-binder propellants, but is 
a must for homogeneous propellants, 

REPEAT FOR NP3 CARDS TO DEFINE ACTIVE BINDER BURNING RATE VS PRESSURE 

LAST CARD 


I 
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B-4 OUTPUT 


The initial portion of the output displays the input as groups labeled 
computation parameters, gas parameters, propellant parameters and the tabular 
input parameters- The ansv/ers are displayed in columns consisting of time, 
web burned, pressure, burn rate, propellant surface temperature, melt layer 
thickness, crater penetration depth, fraction of exposed nitramxne surface area, 
burn area and effective burn rate. Except for bvirn area, the last five parameters 
are associated v;ith nitramine propellants only. The effective burn rate is the 
actual burn rate multiplied by the ratio of the cratered surface burn area to 
the ideal geometric burn area, and would be the value of burn rate deduced from 
a closed vessel test in the absence of knowledge of the actual (cratered) burn 
area. The burn area which is output is the geometric (uncratered) burn area. 

Answers are displayed in time increments of 0,1 msec* 

Computer run time may be controlled by a time stop or max time in 
accordance with particular facility usage, and is usually included in a control card 
The run will automatically stop when the web burns o\it or the pressure of 40000 psi 
is exceeded. This maximum pressure appears in subroutine PRINT, and could be 
changed if desired. Full runs have ranged from 5 min. -15 min* CPU on a 
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